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PREFACE 



As part of the Ministry of the Environment's contribution 
to the International Hydrological Decade, the River Basin 
Research Section is carrying out studies of water resources 
and physical conditions in five basins in Southern Ontario. 
The basins were selected as being representative of common 
type areas in the province and the hydrologic studies being 
undertaken are designed to provide a better understanding of 
most aspects of the water balance in these areas. 

This paper outlines a proposed approach towards the 
development of mathematical models that will describe the 
complex interactions of the various components of the hydrologic 
cycle in each basin. Such models will lead to a better under- 
standing of the actual hydrologic processes that exist and will 
be useful in the synthesis of hydrologic characteristics in 
ungauged areas thereby aiding in the solution of water resources 
problems . 

It is anticipated that the compilation of pertinent 
information and approach outlined in this paper may be helpful 
to others engaged in the development of mathematical hydrologic 
models. 




K. E. Symons , Director, 

Water Quantity Management Branch, 



Toronto, November 1, 1972. 
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AN APPROACH TO MATHEMATICAL MODELLING OF 
MINISTRY OF THE ENVIRONMENT 
IHD REPRESENTATIVE BASINS 

INTRODUCTION 

According to the 'Guidelines for Research Basin Studies', 
as compiled in 1966 by the Canadian National! Committee for the 
International Hydrological Decade, 'representative basins are 
areas set aside to permit study of the existing conditions of 
a region'. One of the objectives of research in these areas 
is: 

'To provide data, on all parameters in 
the hydrologic cycle, so that statistical 
or physical models may be developed for 
major climatic and geologic areas, * 
Principles derived from these models will 
be used to synthesize the hydrological 
characteristics of ungauged areas. 1 

An aim of hydrologic research is to interpret collected 
data for use in the solution of water management problems. 
All drainage basin studies should therefore be planned to 
further the understanding of the hydrologic cycle by using 
knowledge of the physical and mathematical relationships 



* See Appendix I for location maps of basins representative 
of major climatic and geologic areas in southern Ontario, 
under study by the River Basin Research Section. 
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between the various components . The development of 'models' 
from these relationships should then provide the ability to 
accurately reconstruct past and present conditions and to predict 
future events, relating to hydrological phenomena,' These models 
can therefore be useful aids to solving water resources problems 
dealing with evaluation, design, management and regulation, 
Hydrologic Models 

Runoff from homogeneous drainage areas and time-variant flows 
in open channels are affected by factors which are generally 
known and can be defined explicitly. These factors are difficult 
to evaluate in a basin-wide study where data are usually scarce 
and vary both in time and space. However, the analysis of a 
rainfall-runoff process in a complex drainage basin can be made 
feasible by the use of models based on concepts of the system. 
Basin models may be classified into three general categories: 

a ) phy s leal mo de 1 s , 

b) analog models, 

c) mathematical models. 

In the selection of a type of model, consideration must be given 
to cost of development, operational problems, accuracy, inter- 
pretation of results and subsequent model flexibility. 

Physical models are used to represent an idealized physical 
process or environment and are designed for a specific purpose, 
at a specific location. Physical models are employed in 
hydrology to simulate such processes as infiltration, ground- 
water flow and rain fall- runoff phenomena, With minor modifications, 
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these models can often be used for extensive periods of time, 
justifying their cost of construction and operation through the 
importance of such economic essentials as water-supply, flood 
control and navigation, 

Analog models are commonly used to a varying degree in 
solving ground-water flow problems. The principle involved in 
this type of modelling utilizes the analogy between the flow of 
electricity in an electronic circuit and the flow of water in 
an aquifer. Analog models, like physical models, are usually 
designed for a specific location, with the costs being justified 
primarily on a singles-purpose basis. 

Mathematical models allow natural hydro logic phenomena to 
be investigated as systems, The systems approach enables the 
development of quantitative relationships between the various 
hydrologic variables which are usually expressed through 
mathematical formulations (models) for simulation purposes. If 
the hydrologic variables are considered to follow physical laws 
rather than laws of probability, the model is described as 
parametric or deterministic. If laws of probability are considered 
in characterising the hydrologic variables, the model is described 
as stochastic or probabilistic. 

Parametric Modelling 

In preliminary hydrologic studies, parametric models of 
surface runoff usually related runoff to rainfall and drainage 
area, and took the form of a purely empirical relationship or 
of a design procedure following the rational method. With the 
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increase of knowledge in the mechanics of physical hydrologic 
processes, modelling techniques become more complex; hence, 
with a more comprehensive approach, parametric modelling can 
be treated on the basis of component modelling and integrated 
system modelling, 

(a) Component Modelling: 

The overall hydrologic cycle can be divided into several 
component processes such as precipitation, snowmelt, soil 
moisture flow, evapo transpi ration , direct surface runoff, 
base flow and total runoff. Each component or combination of 
components, may be treated as a separate model, For each model, 
empirical approximations have been developed to simulate the 
processes controlling the component (s) , However, by gaining 
an improved understanding of the physical laws governing a 
component, improvements on the initial empirical models have 
been achieved whereby the governing concepts underlying the 
physical situations can be more closely approximated. 

(b) Integrated System Modelling: 

Through the development of better component models, the 
overall modelling of the integrated or total system of a watershed 
has become more complex. Through the use of electronic and 
digital techniques, the limitations on the complexity are 
alleviated and the individual component models can be linked 
to form a comprehensive, integrated system model. 
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Stochastic Modelling 

The approach in stochastic modelling is to present a 
statistical simulation of a measured response of a sy stem- 
In order to use effectively a stochastic model, a probability 
distribution for the input must be determined or assumed. In 
past hydrologic studies, considerable effort was made to 
determine the probabilistic properties of hydrologic data, 
particularly those of streamflow and rainfall. 

One of the characteristic features of the stochastic 
method is that the statistical measures of the hydrologic 
variables used in the calibration of the model enable probability 
limits of accuracy to be placed on the simulated values. 
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OBJECTIVES 

The general purpose of this paper is to present the 
conceptualization in hydrologic modelling relative to present 
and projected applications to hydrologic investigations in 
Ministry of the Environment - IHD representative basins. The 
schematizations developed will serve as flexible frameworks 
for a subsequent detailed approach and a step-by-step procedure 
in the de ve 1 op me n t of a pp 1 i c ab le h y d rol og i c mode 1 s , 

The encompassing objective of the hydrologic simulation 
studies is to formulate models which possess detailed and 
specific features which are characteristic of a representative 
study area. The models that are formulated should possess 
equally sufficient generality for extrapolation to other study 
areas or drainage basins within the larger study region - 
Detailed information, hydrologic and geologic data, compiled 
for the selected representative basins, are to be used in the 
optimization of the model parameters and the overall model 
calib ration. 
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HYDR0L0G1C MODELLING 

The review of methods, approaches and procedures in hydro- 
logic modelling, which are to be used as bases and guidelines 
in the development of working models for each of the Ministry 
of the Environment - IHD representative drainage basins, is 
summarized in the appendices, These summaries are outlined in 
a series of figures as model schematizations . The following 
text serves to introduce briefly the content of each Appendix. 
Accompanying each Appendix is a further elaboration on the text 
which serves to explain in more detail the schemes of the 
respective figures. 

These hydrologie modelling procedures are outlined on the 
basis of: 

1. Basic Conceptualization 
exemplified by, 

i) A schematization for functional and inter- related 
hydrologie processes. 

ii) A schematization for the superposition of basin 
instrumentation on the functional and inter- 
related hydrologie processes. 

2. Component Models 
exemplified by, 

i) A schematization for hydrologie component models 
for data generation and simulation. 
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3. An Integrated Hydrologic System Model 
exemplified by, 

i) A schematization for an integrated (comprehensive) 
hydrologic system model, 

4. A Specific Model Approach 
exeropl i fie d by , 

i) A definition scheme for a parametric, linear 

regression model, 
ii) Expanded definition scheme for a parametric, 

linear regression model, 
iii) Flow chart for a parametric, linear regression 

model . 

5. A Brief on Stochastic Models 
exemplified by, 

i) A schematization for a type of stochastic 
hydrologic system model. 

Basic Conceptualization 

The basic conventional concepts of the hydrologic cycle 
were used to formulate the schematization of Figure 1 of 
Appendix II. The functional and inter- related hydrologic 
processes are outlined by means of a flow scheme showing 
basically the paths precipitation, as input, would follow 
through the schematization to the various outputs. 

The standard, measured meteorological parameters, precip- 
itation, P, temperature, T, radiation, R, — (longwave radiation, 
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R, , shortwave radiation, R , or net radiation, R ) , duration 

J. S il 

of sunshine, S, humidity, H, wind, W, and Class-A pan eva- 
poration, E , are shown as linkages between the source and 
the processes for which they are a cause or effect, or on 
which they have an influence directly or indirectly. 

The precipitation process resulting in rain and/or snow 
en route to supply the surface storages, Liquid Amount and 
Snowpack, is affected by the physical process of Interception 
which results in a reduction in the total amount available to 
these surface storages. In the process, an amount is trapped 
by the vegetal canopy as Interception Storage. The surface 
storage, Snowpack (snow), is subsequently affected by certain 
physical interacting processes (heat exchange and melting) , 
resulting in equivalent liquid amounts. The momentary storage 
at the surface as rainfall amounts and the melt drained from 
the snowpack are the initializing sources of input to the 
Soil Profile Storage, a subsurface storage, which is affected 
by the physical process of Infiltration, 

That portion of the initializing input which is in excess 
of the infiltration requirement, goes to initialize the surface 
physical process, Surface Runoff (overland flow) by way of 
satisfying the Depression and Detention Storages, Eventually, 
amounts from the Depression and Detention Storages, infiltrate 
to the Soil Profile Storage. 
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A portion of the supply to the soil profile, normally 
when in excess of its storage capacity, serves as the source 
for positive supply to the Aquifer Storage, This contribution 
is affected by the subsurface physical process of Percolation, 
which is characterized by a saturated flow. In situations where, 
and at times when, gravity flow from the soil profile to the 
water table does not exist, a reverse flow {unsaturated flow) of 
moisture may be effected by Capillarity; the process being initiated 
and maintained by a negative gradient set up by the effect of 
evaporation and/or transpiration at the soil surface. 

The Aquifer Storage serves as a source for the active ground- 
water flow process resulting in Baseflow into the stream. In 
situations where a local aquifer is in continuum with the expanse 
of a regional aquifer, a portion of the input to the system may 
be accounted for as a loss to deep percolation in the regional 
ground -water flow system. 

Total Runoff, the summation of the Surface Runoff and Baseflow, 
routed by way of channel (s) flow is reproduced at the drainage 
outlet as a Total Runoff Hydrograph. 

Direct losses from the various storages, such as Interception, 
Depression and Detention, Liquid Amount and Snowpack, are affected 
by the process of Evaporation. The Transpiration process by 
plants result in indirect moisture losses from the unsaturated and/ 
or saturated zones . The Evaporation and Transpiration processes 
form, the overall Evapo-transpiration process and the sum of these 
losses constitute a negative output from the hydrologic system as 
Actual Evapo-transpiration, 
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Figure 2 of Appendix II, relates to the various instru- 
mentation and/or methods used for quantitative determination 
of the various hydrologic variables or parameters in terms of 
amounts, durations, intensities and frequencies of measurements. 
For example, the input variable, Precipitation, is sampled 
for amount, duration, intensity and a real distribution. The 
instruments used for the sampling being storage gauges and 
recording and standard rain gauges. 

The scheme for each segment in Figure 2 of Appendix II, 
outlines, in general, the standard instruments used for 
measurement, the units of measurement and the frequencies of 
measurement. In cases where there are no standard instruments 
for use in determining the variable or parameter, the most 
widely used methods for carrying out these determinations are 
indicated. For example, surface runoff may be estimated from 
a total runoff hydro graph by applying one of several hydrograph 
separation techniques. 

Component Models 

Figure 3 of Appendix III shows a schematization for component 
modelling, relating briefly some of the types of models which 
may be investigated under each component, A 'model', in this 
context, is used very broadly; that is, it incorporates all 
operations that result in a historical or synthetic data 
sequence of the variable. The text on component models in 
Appendix III consists of a caption for each model and expands 
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the summarized scheme of Figure 3 f quoting examples of references 
for each type of model, typical data requirements and brief 
comments on present and projected applications relative to the 
River Basin Research Section. 

The block (Figure 3, Appendix III) for 'Data Management' 
outlines briefly the necessary and prerequisite operations to 
be done on the data prior to the development of a component model. 

A group of component models may be combined, through a 
system of transformation {'Transformation System'), to develop 
an overall comprehensive catchment model. The generated or 
historical data sequence for each of the component models may 
form part of the data pool for use in tests for the sensitivity 
and performance of the overall model developed. 

An Integrated System Model 

Figure 4 of Appendix IV, shows a schematization for an 
integrated {comprehensive) hydrologic system model. Outlined 
are a series of operations through various schemes involving 
the Precipitation input as rain, or excess snowmelt, which is 
routed through the various catchment storages, namely, 
Interception Storage, IS, Surface Storage, SS, Aeration Zone 
Storage, AZS, Interflow {intermediate) Zone Storage, IZS, and 
the Saturation Zone Storage, SZS. Routed out of three of the 
storages are outflow components, namely, Overland Flow, OLF, 
Interflow, INFL, and Base flow, BFL. These are subsequently 
routed into Streamflow through a linear or non-linear reservoir 
and channel storage or flow translation procedure to produce an 



- 13 - 

outflow hydrograph Total Runoff, TRQ, at the basin outlet, 
Some of the basic mathematical expressions that relate the 
various parameters within a given storage function are cited 
for the intermediary schemes of operations between the input 
and output. Several explicit relationships have been developed 
by several investigators and were cited by them to be operational 
in an integrated system model through empiricisms and numerical 
approximations. It should be noted, however, that the parameters 
defined in the various functions have to be determined and 
optimized on the basis of sufficient physical justifications. 
The attached text to Figure 4 of Appendix IV, serves to 
define the various parameters of the given functions and outlines 
further a brief account of the operational sequences throughout 
the system. It is not the intention at this stage to elaborate 
on the utilization of these functional relationships, in their 
explicit forms, in actual numerical or computational procedures, 

A Specific Model Approach 

Figures 5a to 5c of Appendix V, are representations of 
basic schemes for a more detailed and specific approach towards 
a framework for a type of parametric model. Starting from a 
simple concept and going through some basic definitions, the 
model is built up to an operating format. Figure 5a of 
Appendix V, is a definition scheme for a parametric model 
(Diskin, 1970) for use in an interpretation of a linear regres- 
sion relationship between precipitation and runoff. Figure 5b 
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of Appendix V, shows an expansion of the model , taking into 
consideration implicitly , retention storage and evapc -transpiration. 
Figure 5c of Appendix V, is the flow chart for the parametric, 
linear regression model outlined in Figure 5a and 5b of Appendix 
V, Given is the flow of mathematical operations and computations 
that may go into the simulation, assuming that the necessary 
parameters, constants and coefficients (regressions and implied 
basin characteristics) were previously determined from other 
regression sub-routines or by other estimating methods. The 
text on the lumped-parametric model of Appendix V, outlines further, 
the step-by-step development of the model, with definitions for 
the: various parameters, coefficients and constants, with an 
explanation of their implied physical significance. Further 
treatment in the text presents a section on the relative stat- 
istical criteria to be used for simplification of the model (e. g, 
operation in one or two elements) for an appropriate application 
in the regression sub-routines. 

A Brief on Stochastic Models 

Figure 6 of Appendix VI, shows a framework for a stochastic 
hydrologic systems model (Chow and Kareliotis, 1970), outlining 
schematically the inter-relationship (formulated on the principle 
of conservation of mass) among the individual models, each of 
which is a component stochastic model. A very long sequence 
of historical record is required for each time series in order 
to adequately determine its probability distribution, or to 
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determine by corzelogram or spectral analysis, the type of 
stochastic model (autocorrelation, harmonic, moving average, 
etc.) to be used in the generating processes. The text of 
Appendix VI gives a short review on an approach to stochastic 
hydro! ogic system modelling. 
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APPLICATION 

The review on conceptual hydro logic modelling outlines 
broadly the extent to which the River Basin Research Section 
plans to proceed with the formulation (building and calibration) 
of a hydrologic model for each of the Ministry of the Environment - 
IHD representative drainage basins. 

The literature reviewed indicates that not many successful 
attempts have been made to formulate one overall hydrologic 
system model which has a general application to all types of 
drainage basins and hydrologic situations. However,, the principles 
and approaches that have been developed from the application of 
modelling to more specific conditions, provide for the selection 
of an approach that may be transposed to areas with similar, basic 
hydrologic characteristics . 

The emphasis in modelling has been on separate hydrologic 
phases (component models). Component modelling, theoretically, 
requires a substantial amount of information and lengthy sequences 
of data for a meaningful approach towards formulating a simulator. 
Also, long sequence of record is required for a totally effective 
calibration of the model. The problem of lack of sufficient data 
and length of records, however, is often alleviated by making 
implicit or explicit assumptions regarding a parameter, group 
of parameters and/or some aspects of the inter- reacting hydrologic 
processes. These assumptions are often desirable for simpli- 
fication of the model to make it computationally tractable. 
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The present extent of the instrumentation in the River 
Basin Research Section's drainage basins and the degree of 
data collection activity, relative to the length of historical 
record of the various hydrologic and geologic variables compiled 
to date, will necessarily dictate the nature and types of 
hydrologic modelling that can be attempted. The relatively 
short sequences of data collected subject the choice of models 
to the confines of parametric and deterministic types and of 
a water balance analysis or accounting approach. Attempts can 
be made as first approximations, at model calibration through 
parameter optimizations, based on averages, ratio estimates or 
least-square methods. 

For the present time, component modelling will be restricted 
to manipulations with the parametric and deterministic types, on 
the basis of tried approaches and existing modelling techniques. 
Knowledge gained from the insight into and/or trials with 
existing model techniques will be used in the eventual formulation 
of an overall comprehensive catchment model, based on the data 
from one or more of the representative basins. 

Although cursory and trial attempts will be made with a 
modelling exercise based on conventional probabilistic treatments, 
no extensive attempts will be made at this stage with stochastic 
modelling due to the lack of extensive data for each component 
time series. Although long series of records are a prerequisite 
for effective utilization of this approach to modelling, attempts 
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may be made to formulate and test simpler existing models, e.g., 
the formulation of a precipitation model generator (Grace & 
Eagleson, 1967; Perkins, 1971), by utilizing relatively short 
sequences of recording precipitation data to develop a stochastic 
model for generating short- time -interval sequences of rainfall 
data. 

Current activities in the Section are centred around data 
management (data reduction, compilation, error analysis, etc) 
which is a prerequisite for a rational and concerted approach 
towards overall water balance analyses and hydrologic modelling. 
Models will initially be of a parametric type, based largely on 
the framework set out in the definition scheme for a parametric, 
linear regression model. Subsequently a more detailed approach 
will be followed, as outlined for an integrated hydrologic 
system model, as shown in Appendix IV. 
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APPENDIX I 



IHD PROJECT REPORTS 



East & Middle Oakville Creeks 
IWB - RB - 16; ONT - 4 

Venison Creek 
IWB - RB - 17? ONT - 4 

Blue Springs Creek 
IWB - RB - 18; ONT - 4 

Bowmanville , Soper & Wilmot Creeks 
IWB - PB - 19; ONT - 4 

Wilton Creek 
IWB - RB - 20; ONT - 4 



Location and Instrumentation Maps of 
Ministry of the Environment- 
IHD Representative Basins 
(in pocket at back of report) 



IHD PROJECT REPORT 



Project Title: Representative Basin Studies - East and Middle 

Oakville Creeks. 
File Reference: IWB - RB ■ - 16; ONT - 4 

Objectives: Collection of hydrologic and geohydrologic data 
necessary for the study of all aspects of the hydrologic regime 
and the general water balance in a basin representative of clay 
and clay till plains in southern Ontario, 

Instrumentation: 5 recording and 3 non- recording streamflow 
gauging stations; 1 main meteorological station, 6 satellite 
meteorological stations; 2 recording, ground-water observation 
wells; 21 peizometers in 8 nests, 12 private, non- recording, 
shallow, ground-water observation wells, 9 snow courses and 12 
soil moisture index stations. 



IHD PROJECT REPORT 

Project Title: Representative Basin Studies - Wilton Creek 

File Reference: IWB - RB - 20 j ONT -4 

Objectives: Collection of hydrologic and gebhydrologic data 

necessary for investigations into all phases of the general 

water balance in a basin representative of limestone conditions 

in southern Ontario. 

Instrumentation: 2 recording and 2 non- recording streamflow 

gauging stations; 2 main meteorological stations and 5 satellite 

meteorological stations; 4 recording, ground-water observation 

wells and 11 piezometers in 4 nests. 



IHD PROJECT REPORT 



Project Title: Representative Basin Studies - Bowman vi 1 le , 

Soper and Wiimot Creeks 
File Reference: IWB - RB - 19; ONT - 4 

Objectives: Collection of hydrologic and geohydrologic data 
necessary for the study of all aspects of the hydrologic 
regime and general water balance in a basin representative of 
moraine, till and clay areas in southern Ontario. 
Instrumentation: 16 recording and 2 non- recording streamflow 
gauging stations; 2 main meteorological stations and 12 satellite 
meteorological stations; 12 recording ground-water observation 
wells; 17 piezometers in 8 nests? 23 private, shallow, non- 
recording ground-water observation wells; 16 soil moisture 
index stations; 2 ground- temperature recorders and 12 snow 
courses . 



IHD PROJECT REPORT 



Project Title: Representative Basin Studies - Venison Creek 

File Reference: IWB - RB - 17; ONT - 4 

Objectives: Collection of hydrologic and geohydrologic data 

necessary for the study of all aspects of the hydrologic regime 

and general water balance in a basin representative of sand 

plain conditions in southern Ontario, 

Instrumentation: 3 recording and 1 non- recording streamflow 

gauging stations; 1 main meteorological station and 3 satellite 

meteorological stations; 4 recording, ground-water observation 

wells and 3 piezometers in 1 nest, 



IHD PROJECT REPORT 



Project Title: Representative Basin Studies - Blue Springs Creek 
File Reference; XWB - RB - 18; ONT - 4 

Objectives: Collection of hydrologic and geohydrologic data 
necessary for the study of all aspects of the hydrologic regime 
and general water balance in a basin representative of kame, 
till and bedrock conditions in southern Ontario. 
Instrumentation: 3 recording and 6 non-recording streamflow 
gauging stations; 5 recording ground-water observation wells ; 
11 piezometers in 3 nests and 8 private, no n- re cording, shallow 
ground-water observation wells. 
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HYDROLOGIC MODEL : 

Basic Conceptualization 

Figure 1 shows the schematization outlining the functional 
and inter-related processes based on the conventional concepts 
of the hydrologic cycle, 

The end product of the meteorological process, Precipitation, 
is rain and/or snow, a portion of which is intercepted by the 
vegetal canopy as Interception Storage and eventually lost to 
Evaporation (or Sublimation) . Precipitation as a liquid is held 
momentarily at the ground surface as Detention Storage. Some, 
through the process of Infiltration, enters the Unsaturated Zone 
and that in excess of the soil moisture capacity goes to replenish 
the local ground-water flow system (Saturated Zone) by Percolation. 
Precipitation with rates in excess of the infiltration capacity 
at the soil surface contributes to Total Runoff by way of direct 
Surface Runoff. A portion of the precipitation to direct Surface 
Runoff is trapped as Depression Storage in ponds, swamps, etc. 
Subsequently, this retained portion of the runoff contributes 
respectively to Evaporation, to Infiltration to the Soil Profile 
and recharge to the local Aquifer. Portions of the recharge 
eventually contribute as Baseflow to Total Runoff and to Regional 
Flow as Deep Percolation, Eventually, there may be a reverse 
flow to the Unsaturated Zone as a result of a capillary (suction) 
potential set up by Evaporation and/or Transpiration at the soil 
surface. 

Precipitation which occurs as snowfall is stored on the 
ground as Snowpack. The pack undergoes various forms of met- 
amorphosis prior to the active snow-melting process. During 
the metamorphosis phase, a quantity of moisture is lost to 
Evaporation (by Sublimation) ? also during this period, several 
phases of heat exchange take place at the snowpack-soil interface 
and between the snowpack and the liquid rain, These and several 
other types of energy fluxes (Radiation, Temperature) influence 
the onset of the Snowmelt process by way of convective heat 
transfer between the pack and atmosphere. The Snowmelt, as 
liquid excess drained from the Snowpack, follows a similar path 
to that of the rain; that is, via Depression and Detention 
Storages and Infiltration, and subsequently to Surface Runoff 
and Baseflow to Total Runoff which are reproduced as a Total 
Runoff Hydrograph at the basin outlet. 

Transpiration by crops from the Unsaturated and/or 
Saturated zones, together with evaporation from the various 
processes (Interception, Snowpack, Detention and Depression 
storages) , contributes to the Actual Evapotranspiration, 
from the basin. 
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The ancilliary meteorological variable Temperature, T, 
is associated with Radiation, R, (as longwave, shortwave or net 
radiation) , with net energy (heat) flux at the incident surface. 
The Humidity, H, is a function of the vapour pressure, e a and 
the saturated vapour pressure, e s , at a given air temperature. 
The Wind, W, is expressed in terms of the wind speed, u, and 
direction. Sunshine duration, S, is directly related to the 
temperature and shortwave radiation variables. These ancilliary 
meteorological variables are necessary prerequisites for use in 
computation (through empirical formulae) of snowmelt and eva- 
poration from snowpack and evaporation from available sources of 
free water. 

Instrumentation 

Figure 2 shows the various instrumentations and/or methods 
used for quantitative determination of the various hydrologic 
variables or parameters. 

The rainfall is measured by standard rain gauges, RS, 
providing total depth quantities in inches, arranged to yield 
daily, monthly and yearly summaries. Rainfall intensities 
and durations are measured by recording rain gauges, RA, providing 
summaries of storm amount, duration and rates per unit time which 
are the data necessary for such analytical procedures as storm 
intensity- frequency and depth-area-duration analyses. Snowfall 
may be measured accumulatively with Alter-shielded or Nipher- 
shielded storage gauges; accumulated snowfall on the ground as 
snowpack is core-sampled for depth, water equivalent and density. 

The other meteorological variables are measured by con- 
ventional instruments: temperature extremes by maximum and 
minimum thermometers; humidity and temperature by hygrothermo- 
graphs; duration of sunshine by Campbell-Stokes solarimeters; 
longwave and shortwave (Solar) radiation by Eppley or Kipp 
pyranometers; net radiation by CSIRO pyrradiometers; evaporation 
by Class-A pans and wind speed by anemometers, 

Interception (amount of precipitation intercepted by the 
vegetation) may be determined by plot studies, or estimated with 
the use of existing, general empirical equations, 

Snowmelts are usually determined by computations from 
snowmelt indices and/or with general empirical equations. In 
the synthesis of snowmelt hydrographs , some basic factors have to 
be determined. These are area! extent of snow cover; snowpack 
parameters with respect to depth, water equivalent, temperature 
and free water; the variation of these parameters with elevation; 
albedo of the snow surface; and the critical sequence of factors 
affecting melt such as melt rate and loss and runoff conditions 
of the basin. 



* 
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The extent of the depression storage capacity and its 
distribution may be determined from topographic maps and field 
surveys of the drainage area. 

The infiltration function may be approximated by empirical 
procedures such as sample plot studies (infiltrometer, rain- 
ulators) , or by laboratory experiments with field soil samples. 

A time series of soil moisture data may be determined by 
thermo gravimetric analyses of field samples or monitored in 
situ by tensiometers or electrical resistance block methods or 
by the nuclear (neutron meter) method, etc. Soil temperature 
profiles may be monitored in situ by thermohms or by ground- 
temperature recorders. The necessary soil moisture characteristics, 
(soil moisture-hydraulic conductivity and soil moisture-tension 
curves) , for describing the unsaturated flow, may be determined 
through laboratory experiments on specific soil samples. 

Aquifer characteristics (hydraulic conductivity, transmis- 
sibility) may be determined by pumping tests and by analytical 
techniques using ground-water level fluctuations in observation 
wells and piezometers. The saturated hydraulic conductivity 
characterizes the rate of percolation. 

Total runoff may be determined by the use of conventional 
streamflow gauging stations (recording gauges with natural 
controls, NR, or artificial controls, AR and staff gauges, NN) , 
Measured stages at each gauge are converted to continuous 
discharge as mean daily flows in cubic feet per second, with 
the use of empirically developed stage-discharge rating curves. 

Base flow and direct surface runoff may be determined by 
analytical methods such as hydrograph separation. A direct 
estimate of base flow from some areas may be determined by the 
measurement of spring discharges. 

Estimates of the potential evapotranspi ration may be 
obtained from adjusted evaporation pan measurements or estimated 
with the use of existing empirical formulae. Estimates of actual 
evapotranspi ration may be obtained from lysimeter studies or 
from adjusted potential evapotranspiration with a crop consumptive- 
use factor determined from land-use studies. 
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Figure 1. HYDROLOGIC MODEL: 



A Schematizatian for" Functional and 
Inter-related Hydro! ogic Processes 



Precipitation 

Is storage gauges, standard 
rain gauges, recording 
9a ugas , etc . 
UM: inches , ins /hour , ins /day 
Ft continuous, hourly, daily, 
year around 



Interception 
- rational estimates 
(empirical equations 
based on plot studies)* 



Rain 

I: standard £ recording rain 

gauges, Rs, RA 
UK: inches, ins /hour, ins /day 
Pi continuous, hourly, daily, 

year around 



Sunshine 

1 1 solarimeter 
UM: hours/day 

Ft continuous, daily, 
year around 



Temperature 

1 1 maxynum-mi nimuni 
thermometers 

UM: °F, °C 

P: daily, year around 



Radiation 
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UMi Langley leal /cm /min.) 

Ft continuous. Integrated 
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Humidity 
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OMi rel, humidity (1) , 
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HYDROLOGIC MODEL ; 

Component Models 

In general, a component model refers to a mathematical 
formulation, involving operations with a data sequence of a 
single hydrologic process. This may be used for extending the 
data sequence or used to produce synthetic quantitative traces 
of the variables for use in simulation. These component models, 
schematized in Figure 3, are treated on the basis of lists of 
different methods and existing procedures that may be used for 
a model development. Quoted are some of the references relating 
to types of models. 

The likely applicability of a type of model (s) to River 
Basin Research (RBR) present and projected investigations in 
the representative basins is also briefly discussed. 

The cited references are given as examples of different 
types of investigations relating directly or indirectly to 
modelling. An attempt was made to identify the reviewed 
literatures and group them under the different component models. 
Identification was based on the type of hydrologic process and 
the choice of the method used in the modelling. 

Precipitation Models 

Method : De te r mini s t i c 

i) inter-station relationship (correlation and 
regression models) for use in data supple- 
mentation and extension, 
ii) storm transposition. 

Probabilistic 

i) order statistics, e.g. conventional procedures 
for depth- are a- duration analysis. 

ii) frequency analyses (probability of discrete 
events) , e.g. storm intensity -frequency 
analysis and extreme values analysis, 

iii) random sequences, (stationary and non- 
stationary process) , e.g. stochastic models 
for synthesis of hourly rainfall data, models 
of the types formulated by Grace and Eagleson 
(1967)*, Pattison (1965) , Franz (1971), 
Raudkivi and Lawgun (19 70) , Hiemstra and 
Creese (1970) . 



* References appended in Bibliography. 
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Data Requirement 

- historical data: daily, monthly and yearly summaries from 
the basin network stations; recording precipitation gauge 
data summaries? monthly, yearly, summaries from selected 
regional index precipitation stations, 

Projected RBR Application 

- each basin has 1 or 2 recording rain gauges and a network 
of standard rain gauges which provide the required data 
summaries . 

- historical records of several long-term regional stations 
within each climatic region are available from monthly 
climatological station reports, published by the Atmospheric 
Environment Service, Canada, 

- application in modelling may be attempted with data extensions, 
order statistics and sample computations based on frequency 
analysis and probability of discrete events. Some attempts 
may be made with stationary stochastic methods. This involves 
utilizing a time series of recording rain gaugedata to 
formulate a stochastic model for use in simulating short- 
sequence rainfall data. 

Snowmelt Models 

Method: Deterministic 

i) empirical equations - snowmelt synthesis 

with the use of general empirical snowmelt 
equations of the Corps of Engineers, U, S, 
Army (1960). 

ii) melt indices - use of melt indices, _ e.g. 

degree-day index, or melt-plot studies, model 
by Watt and Hsu (1971), 

iii) index plot and empirical equations - 

snowmelt simulation using a hybrid of index 
plot studies and empirical formulae, e. g, 
models formulated by Anderson and Crawford 
(1964), Amorocho and Espildora (1966), 
Crawford and Linsley (1966) . 

Data Requirement 

- basic meteorological parameters (rainfall, temperature, 
radiation, vapour pressure, duration of sunshine, wind, 
convective heat transfer) , 

- areal extent of snow cover and snowpack parameters (depth, 
water equivalent, density, temperature, albedo of snow 
surface, permeability, etc.). 

- streamflow records. 
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Projected RBR Application 

- basic meteorological parameters are monitored continuously 
in each basin. 

- a real extent of snow cover and snowpack parameters (depth, 
water equivalent, density, etc.), are determined seasonally 
in two basins, 

- winter streamflows are monitored continuously. 

- present application in modelling may be attempted with 
snowmelt hydrograph synthesis utilizing general melt 
equations and/or melt indices, e. g. models of the type 
formulated by Watt and Hsu (19 70) which utilize a melt- 
index based on the degree -day concept, or a modified approach 
to the models used by Amorocho and Espildora (1966) which 
utilize empirical equations based on energy-melt concepts 
and index plot observations. 

Soil Moisture Models 

Method : Deterministic 

i) * lumped- system as related to conceptual 

catchment storage; the use of soil moisture 
data, for specified soil types, in a budget- 
type equation to develop empirical relation- 
ships for use in estimating soil moisture 
regimes , e.g. models of the type used by 
Zahner (1967), Visser (1966), Baier (1969), 
Baier and Robertson (19 65) . 
- the use of depletion curves and accretion 
regressions from measured field data to 
develop empirical relationships for predicting 
moisture regimes, e. g. models of the types 
used by Carlson et al (1956) , Stearn and 
Carlson (19 60) . 

ii) lumped-system as related to infiltration; the 
use of soil moisture parameters in the 
modelling of infiltration process with methods 
such as those used by Overton (1964), 
Muggins and Monke (1968) , Deboer and Johnson 
(1969), Skaggs et al (1969), Ligon, et al (1965). 

Data Requirement 

- soil types, land use and cropping practices and their areal 
distribution. 

- hydrologic characteristics (moisture retention curves, 

* 'lumped-system' - the spatial distributions of the input and 

output are ignored. 
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infiltration rates, hydraulic conductivity, etc.) of the 
major soil types. 

- historical data of absolute soil moisture or changes in 
soil moisture for the study area, point measurements from 
a network of stations. 

- meteorological variables (rainfall, temperature, radiation, 
humidity, evaporation, etc.). 

Projected RBR Application 

» soil moisture observations (absolute and changes in moisture 
from field studies ) are being carried out in one sub-basin 
(BSW-Wilmot) and the initiation of a similar study is in 
progress for EMO . 

- summaries of major soil types and land use and cropping 
practice in these basins are compiled. 

- soil hydrologic characteristics to be determined by lab- 
oratory experiments on field samples. 

- present application to modelling may be attempted with 
deterministic methods, such as the use of a water balance 
type equation to develop empirical relationships from field 
observations for use in estimating soil moisture regimes, 
e.g. method by Zahner (1967) , Baier and Robertson (1965) ; 
the use of correlation relationships between soil moisture 
changes and meteorological variables to predict moisture 
regimes, e.g. method by Steam and Carlson (I960); and the 
use of soil moisture parameters in infiltration models, e.g. 
method of Huggins and Monke (1968). 

Evapo transpiration Models 

Method : Deterministic 

i) empirical relationships for potential evapo- 

transpi ration; the use of developed evaporation 
formulae which are based on the Energy-Budget 
and Mass-Transfer theories of evaporation such 
as the Penman or Thorn thwaite formulae and 
method by Gay (1970); the use of adjusted 
Class-A pan data with other meteorological 
variables to develop correlation or regression 
relationships for use in predicting potential 
evaporation regimes, such as done by Mukammal 
and Bruce (1960), Christiansen (1966), Morton 
(1971) , Harmon (1966) . 
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ii) empirical relationships for determining actual 
evapo transpiration; the use of a water balance 
type equation with measured field data to 
develop an empirical relationship for use in 
predicting the evapotranspi ration regime, 
such as the method used by Bowman and King 
(1965) ; the use of a meteorological budgeting 
technique with meteorological and soil moisture 
data to develop empirical relationships for 
estimating daily evapo transpiration, e.g. 
method by Baier (1966) % the use of pan eva- 
poration and meteorological variables to 
develop relationships for estimating actual 
evapotranspi ration such as summarized by 
Hargreaves (19 60) ; the use of solar radiation 
and temperature to estimate evapo transpiration, 
e.g. method by Robb (1966). 

P robab il i s t i c 
i) time-series model? the use of time-series 

analysis on historical sequences of estimated 
potential evapo transpiration to determine a 
stochastic model for data generation and 
simulation, Chow and Kareliotis (1970) , 

Data Requirement 

- historical data: hourly, daily, monthly or yearly summaries 
of meteorological data (radiation, temperature, humidity, 
wind speed, duration of sunshine, Class-A pan data, etc.). 

Projected RBR Application 

- the basic meteorological variables including Class-A pan 
evaporation data are monitored in each basin, 

- application to modelling may be attempted with computations 
based on available empirical formulae or by trial computation 
with formulations through correlation and regression analysis, 
e. g. the application of a modified Penman's or Thornthwaite 1 s 
equation and the development of regression relationships 
between Class-A pan evaporation and meteorological variables 
for use in data supplementing or data extension, or the use 
of a budget equation with measured meteorological and soil 
moisture data to develop relationships for use in estimating 
actual evapo transpiration. 

Surface Runoff Models 

Method ; Deterministic 

i) linear time-invariant; the use of a distributed 
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linear method for the synthesis of runoff 
hydrographs from given rainfall inputs, such 
as the Unit Hydrograph Method, as applied by 
Eagleson (1962), Mitchell (1967), Gray (1961), 
Viessman (1968), Overton (1967), Dooge (1959) ? 
Eagleson et al (1965) ; or the use of other 
routing methods such as the rational method 
and the summing of hydrographs from a 
number of sub- are as within a basin. 

- the use of the physiographic parameters of 
the basin and surface runoff characteristics 
to develop regression equations for use in 
simulation, e. g. as developed by Narayana 
and Riley (1971), Schulz et al (1971), 
Taylor and Schwarz (1952), Larson (1965). 

- the use of a hydraulic model of the catchment 
to simulate runoff, e.g. based on Kenematic 
wave theory as applied by Wooding (1965, 
1966). 

ii) non-linear time -in variant,* the use of a 

distributed non-linear method to synthetise 
runoff hydrographs, e.g. IUH as applied by 
Singh (1964) . 

Probabilistic 

i) frequency analyses; the use of peak flow 
data for annual or partial duration series 
to develop frequency curves and statistics 
for use in flood flow predictions, e. g. 
as applied by Snyder (1958) . 

- the use of the statistical parameters of 
flood series and rainfall extremes to develop 
regression relationships for use in simulation, 
e. g. as applied by Reich (1970), Shane and 
Gaver (1970). 

ii) random sequences (stationary and non- 
stationary processes) ; the use of the 
stochastic property of flood flow series to 
develop a model for use in simulation, e.g. 
as applied by Todorovic and Zelenhasic 
(1970), Shane and Lynn (1964); the application 
of a stochastic model to storm runoff and 
rainfall data to formulate relationships 
for use in generating sequential data, 
after Chow and Ramasseshan (1965). 

Data Requirement 

historical data: streamflow (mean daily and peak flows) ; 
recording precipitation summaries (amount, intensity, duration) . 
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Projected RBR Application 

- the use of deterministic and probabilistic methods - re 

synthesis of surface runoff hydrograph, with the use of unit 
hydrograph or instantaneous unit hydrograph techniques; re 
development of relationships between runoff hydrograph 
parameters and basin physiographic parameters; re development 
of frequency curves and parameter relationships for flood 
flow series } re development of regression relationships between 
flood flow and rainfall extreme series, etc. 

Ground -Water (Base flow) Models 

Method : Deterministic 

i) lumped- system; to consider the relationship 
between total quantity of precipitation and 
its direct runoff component, for use in 
estimating ground-water flow and rainfall 
loss relative to the infiltration capacity 
and soil moisture as applied by Kadoya (1967) , 
Meyboom (1961) . 

ii) distributed- system; to consider the theory of 

flow problems re models based on the analytical 
and numerical solution of the flow equation with 
given boundary conditions for specified problems, 
e. g. for a simplified, two-dimensional flow 
problem as applied by Toth (1962, 1963) and 
for multi-layered , two-dimensional flow 
problem as applied by Freeze and Witherspoon 
(1967); to consider heterogeneity of basin 
models by determining the statistical variability 
of input data (conductivity, specific yield, 
water-level elevation, spatial dimension of 
aquifer, etc.) e. g. solution of the flow 
equation with variable coefficients as applied 
by McMillan (1966), Shahbazi and Todd (1967), 
Kiraly (19 71) . 

Probabilistic 

i) correlation; to develop statistical relation- 
ships between base flow data (determined by 
hydrograph separation techniques) and drainage 
basin parameters or meteorological parameters, 
e. g. relationship between drainage density 
and base flow as applied by Trainer (1969) , 
between base flow and basin parameters as applied 
by Furness and Busby (1969) , between ground- 
water chemistry and base flow as applied by 

* 'distributed-system' - the spatial distribution of the input 

and output is considered. 
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Pinder and Jones (1969) ; relationships of 
quantitative geomorphology to base flow as 
applied by Farvolden (1963) ; relationship 
of recharge to ground water from storm inputs 
to increment in base flow as applied by 
Trainer (1969) . 

Data Requirement 

- surficial and subsurface geology. 

- summaries of water-table fluctuations , spatial dimension of 
aquifers, chemical properties of ground water, pumping tests, 
etc. 

- aquifer characteristics and their variability (specific yield, 
hydraulic conductivity, etc.) 

- natural spring discharge measurements (where possible) . 

- low- flow data (streamflow hydrograph) . 

Projected RBR Application 

- compilation of basic data (as above) is in progress. 

- the use of parametric methods such as the lumped-system to 
relate indices of gross precipitation inputs, direct runoff 
components and corresponding antecedent moisture conditions. 

- the use of simplified one or two-dimensional flow problem 
models . 

- the use of correlation relationships between base flow components 
and drainage basin or meteorological parameters. 

Total Runoff Models 

Method : Deterministic 

i) linear time-invariant; the use of an 

antecedent-precipitation-index type of 
rain fall- runoff relationship and a 
retention index to simulate continuous 
streamflow hydrographs from rainfall inputs, 
e.g. as applied by Sittner et al (1969). 

Probabilistic 

i) regression and correlation; the use of 

multiple regression techniques to develop 
relationships between streamflow data at 
inter- related stations; to develop relation- 
ships between streamflow parameters and 
basin physiographic parameters; to develop 
relationships between the original variables 
or statistical parameters of rainfall and 
streamflow series, as applied by Beard (1965), 
Benson and Matalas (1967), Caffey (1965). 
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ii) auto-correlation; the use of stochastic 

models (Markov Models) on historical stream- 
flow data to develop relationships for use in 
generating synthetic data; the use of a 2nd- 
order Markov model, applied to standardized 
variates (assumed log- Pearson Type III 
distribution), e. g. method by Beard (1967); 
the use of Markov model for low-flow analysis, 
e. g. method by Fiering (1964) ; the use of a 
Ist-order Markov model and regression analysis 
on multiple -station flow data to develop a 
simulator for generating daily flows, after 
Payne and Newman (1969) . 

iii) other time-series analyses; the use of time 

series analyses on historical streamflow data 
to determine the inherent stochastic model for 
use in data generating; the autocorrelation 
and spectral analysis of streamflow series, 
after Quimpo (196 7) , Roesner and Yevdjevich 
(1966), Svandize (1967). 

Data Requirement 

- historical data: streamflow (mean daily, monthly and yearly 
summaries) for both base and secondary gauging stations; 
precipitation summaries* 

Projected RBR Application 

- compilation of streamflow data for each basin gauging network 
is in progress , 

- application of deterministic and probabilistic methods - re 
synthesis of total runoff hydrographs , using a developed 
relationship between antecedent-precipitation-index, retention- 
storage- index and runoff, or using a multi-station regression 
relationship, or using a regression between streamflow and 
basin parameters and/or meteorological parameters, 

- application of probabilistic methods - re use of simpler Markov 
models for data generation (generating models based on mean, 
standard deviation and serial correlation of a flow sequence) . 

Transformation System (Catchment) Models 

Method ; Deterministic 

i) lumped-system; to consider the catchment as a 
linear system, governed by a differential 
equation, the integration of which provides 
the convolution equation for the linear 
transformation of an input (index of precip- 
itation) to an output, as applied by Amorocho 
and Orlob (1961), O'Connell and Nash (1970). 
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ii) distributed-system; to consider the catchment 
to be composed of an infinite array of lumped, 
independent linear subsystems with each sub- 
system being governed by a linear differential 
equation, the summation and integration of 
which leads to the convolution relationship 
for transforming inputs to outputs, as applied 
by Eagleson (1967), Wooding (1965, 1966), 
Brakensiek and Ons tad (196 8) , Muggins and 
Monke (1970), Foster, et al (1968), 

iii) linear reservoir; to consider the catchment as 
a single reservoir or a cascade of n-linear 
re se rvo irs, in e ach of wh i ch the s to rag e is 
linearly related to the outflow, e. g. use of the 
storage equation and the hydro logic continuity 
equation to derive the differential equation 
that governs the linear reservoir, as applied 
by Nash (1957) , Sarma et ai (1969) , Overton 
(1967), O'Connell, et al (1970). 

iv) linear channel and linear reservoir; to 

consider the catchment to consist of a series 
of linear channels and linear reservoirs 
having a combined linear translation effect, 
using a time- area-concentration diagram as 
applied by Dooge (1967) . 

v) non-linear channel and reservoir; to consider 

the catchment storage to be distributed and 
non-linear by dividing the catchment into 
sub-areas, considering their individual 
inputs and routing the excess from each sub- 
area through different amounts of storage, as 
applied by Laurenson (19 64) , Deboer and Johnson 
(196 9) ; to consider the operation of the 
catchment by a functional series (the gen- 
eralization of the U.K. concept to that of a 
non-linear response function) as treated by 
Amo ro ch o and Or lob (1961), M a chme i e r and 
Larson (1967), Overton (1967), Prasad (1971), 

Probabilistic 

i) the use of a linear regression relationship 
between inputs and outputs of the system, as 
applied by Diskin (19 70) ■ 

ii) to consider the catchment as a stochastic hydro- 
log ic system; formulate a system model on the 
basis of conservation of mass regarding the 
components (precipitation, runoff, storage 
and evapotranspi ration) as individual stochastic 
processes, as applied by Chow and Kareliotis 
(19 70) . 
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Data Requirement 

physiographic parameters of the drainage basin (area, 
drainage density, slope, main stream length, etc.). 
historical data or hydrologic parameters (precipitation, runoff, 
etc. ) 

RBR Projected Application 

the use of the lumped-systems re application of a convolution 
integral to transform lumped inputs to outputs. 
the use of simpler distributed system re application of the 
convolution integral to sub-systems. 

the use of a single or double linear reservoir storage re 
application of the storage and hydrologic continuity equations 
to derive the differential equation that governs the input- 
output relationship. 

the use of linear channel translation and linear reservoir 
storage re application of a time-area-concentration relationship 
for the transformation of inputs to output. 

the use of a linear regression relationship between precipitation 
and runoff. 
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HYDR0L0G1C MODEL : 

An Integrated Hydrologic System Model 

Figure 4 shows a schematization for an integrated hydrologic 
system model. 

Precipitation 

Precipitation, P, as rain and/or snow (snowmelt) is the 
initializing input to the system. The initial amount is inter- 
cepted to some degree by the vegetal canopy. 

Interception Storage 

The intercepted precipitation is held by the vegetal canopy 
as Interception Storage, IS. This E is a function of the amount 
of P, the duration of the precipitation, t f the evaporation 

rate, E , at the intercepting surface, the interception storage 
a 

capacity, ISC, of the intercepting area and the ratio of the 

evaporating surfaces to the total horizonally projected canopy 

surface area, A. An explicit function relating these parameters 

is required for use in determining IS from the gross P. 

Surface Storage 

The excess precipitation goes to satisfy the Surface Storage, 
SS, requirement. The Surface Storage is composed of three 
conjunctional functions, namely the Depression Storage, DS , the 
Detention Storage, D, and Infiltration, I. Some explicit function 
is required to distribute the excess precipitation to the require- 
ments of DS, D and I. Generally I reacts simultaneously with DS 
and D to the precipitation input, but DS is normally satisfied or 
partially satisfied before the initiation of D, 

Depression Storage 

The Depression Storage, DS , is a function of the depression 
storage capacity, DSC, of the drainage area, the amount of 
precipitation excess, P , the infiltration characteristics of the 

soil, ic, and the ratio of the depressional storage area to the 
total drainage area, B. 

Detention Storage 

The Detention Storage, D, is a momentary storage which is a 
function of the overland flow length, L, the infiltration 
characteristics, ic, and the overland flow rate, q. The flow 
rate q, in turn, is a function of the overland flow length, L f 
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the slope, S, and the flow characteristics (Manning's roughness 
coefficient, n) of the overland flow path and the rate of 
precipitation input, i. Whenever the input rate, i, is greater 
than the rate of infiltration, overland flow, OLF, is initiated 
and is routed out of Detention Storage D in terms of q, 

•The OLF out from D is routed to the streamflow according 
to some method of linear or non-linear channel storage and flow 
translation to produce the direct Surface Runoff, SRO, component 
of the Total Runoff, TRO, measured at the basin outlet. 

Infiltration 

The Infiltration rate, I, is a function of the soil moisture, 
SM, the amount of excess precipitation, P @ , the rate of precip- 
itation input, i, the initial infiltration rate, I , and the 
infiltration capacity, I , of the soil. Infiltration proceeds 

at a decreasing rate with time; the initial rate is dependent 

on the antecedent soil moisture condition and continues at a 

rate subject to I , until the liquid amounts in DS and D become 
•* c 

limiting. Several empirical infiltration formulae have been deve- 
loped by many investigators (Horton, Hoi tan, Huggins and Monke, 
Phillip, et al) . A few of these empirical formulae have been 
suitably applied to describe the infiltration process in source 
accounting type hydrologic models. The infiltration rate function 
may be examined by the standard flow equation which relates the 
inflow rate, V , to the hydraulic conductivity, K, and the dif- 

fusivity, (D = K/ d (SM ) ,) , of the soil and to the differential 

dip 
change in soil moisture with depth, 3 (SM) , relative to the soil 

3 y 

surface. An empirical determination of K and D and a numerical 
approximation technique for solution, are required for incor- 
poration of the infiltration rate function into a working model. 

Aeration Zone Storage 

Infiltrated water moves through the soil profile in^ 

accordance with known physics of flow through porous media 

thereby satisfying the moisture requirement of the Aeration 

Zone Storage, AZS. The soil moisture, SM, is a function of the 

flow characteristics of the soil, i.e. unsaturated hydraulic 

conductivity, K , capillary potential,^, gravity potential, 

unsEu . 

y and total porosity, T , of the medium. Flow of water in the 

unsaturated zone can be described by the expanded partial 
differential Darcy flow equation, which is the same standard 
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flow equation for infiltration above , but written in the more 

classical form, 3 (SM) m K o . ( If + 1) , relating the rate 

11 unsat. j*~ 

of change of moisture content with the hydraulic potential. 

Interflow Zone Storage 

Infiltrated water, drained by gravity from the Aeration 
Zone Storage, AZS, percolates to the Interflow Zone Storage, IZS. 
Saturated flow is characteristic of the IZS zone. This zone is 
in fact a transient saturated zone or a saturated zone created 
by a less permeable layer below the more permeable soil profile, 
thereby forming a 'perched water table 1 . The physics of 
saturated flow can be used to describe the seepage to be routed 
out of this IZS, as interflow, INFL, which is subsequently routed 
into the streamflow as part contribution to Total Runoff, TRO, 

Saturated Zone Storage 

Drainage by percolation from the IZS goes to recharge the 

Saturated Zone Storage, SZS, Moisture flow in this zone is 

characteristically saturated flow. The outflow q, is dependent 

on the amount Q in storage and the total hydraulic potential, 
o 

H; the outflow rate being regulated by the saturated hydraulic 

conductivity, K . . For saturated flow, SM and K are 

s at • s at • 

constants? and 9{SM)/9t * o. With these known conditions the 
differential flow equation may be solved for unconfined steady- 
state flow in an isotropic porous medium resulting with the Laplace's 
equation ?{KVH) = 0. For a given drainage problem with defined 
boundary conditions, the differential flow equation for the 
problem may be expressed explicitly. With the use of numerical 
approximation techniques , the explicit flow equation may be 
solved approximately and so facilitates the routing of a base- 
flow component BFL, out of SZS, into streamflow as another part 
contribution to TRO. 

The sum of INFL and BFL constitute the subsurface runoff, 
SSRO, component to TRO. 

Total Runoff 

The Total Runoff, TRO, may be routed through the channel 
reach by a linear (or non-linear) storage equation utilizing 
the law of continuity in the basic routing procedure. This TRO 
may be reproduced at the basin outlet as the streamflow hydro - 
g raph . 
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E vapo tr ans p i r a t io n 

Actual Eva po transpiration, ET f is directly proportional 
to the amount of moisture available in the various storages, 
SSM, to the amount of available energy fluxes in terms of 
radiation, R, and temperature, T, to a crop use factor C^ and 

wind speed u, and is indirectly proportional to the relative 
humidity, H, of the environment. 

Evaporation takes place, from the IS and SS at the potential 
rate until the free water available in the sources is used up. 
Evaporation also takes place from AZS, IZS and/or SZS, at a rate 
which is proportional to the rate at which the water is supplied 
to the soil surface through capillarity. 

Transpiration, by plants, from the AZS, IZS and/or SZS is 
controlled by the amount of moisture available to the plant. 
In the case of the AZS, this availability of moisture is limited 
to the range between that of field capacity and wilting point. 

An explicit function (not shown diagr annua tically) such as 
one relating ET to the potential evapotranspi ration rate and soil 
moisture, SM , in the AZS, to crop-use, is required for incor- 
poration into the model to account for ET loss from the basin. 

Working Model 

The development of a working model based on, or modified 
from, Figure 4 requires that a choice be made of an existing 
explicit function for each component of the model. Explicit 
functional relationships may also be developed independently 
through empiricism based on observations and plot studies. A 
linkage or step relationship is also required for the computational 
sequences between the different component schemes. 



LIST OF SYMBOLS IN MODEL 



IS = INTERCEPTION STORAGE - in, 

P = gross precipitation - in. 

t - duration of storm - in/hr. 

E r - rate of evaporation from the intercepting 

a surface - in/hr, 

A ■ ratio of evaporating surface to projectional 

area - sq.mi. 
ISC = storage capacity of the intercepting surface 
areas - in . 

SS = SURFACE STORAGE - in. 

PS ■ Depression Storage - in. 

DSC - Depression Storage Capacity - in. 

ic - refers to the infiltration characteristic 

of the soils 
P ■ excess precipitation (P - IS) - in. 
B e = ratio of depressional storage area to total 

drainage area. 

D m Detention {overland flow) Storage - in, 

L = length of the overland flow reach - ft. 

ic • infiltration characteristic of the soils 

q = rate of overland flow per unit area - in/hr, 

S - slope of the overland flow surface - ft/ft. 

n - a Manning's roughness factor for overland flow 

i - the inflow (precipitation) rate - in/hr, 

OLF = overland flow - in. 

I = Infiltration Rate - in/hr, 



SM m soil moisture -in. 

I ■ initial infiltration rate - in/hr. 
o 

I ■ infiltration capacity - in/hr. 

c 

P ■ amount of precipitation excess - in. 

e 

i = inflow (precipitation) rate - in/hr. 

- Flow Equation - 

V - inflow rate - in/hr. 
o 

K ■ hydraulic conductivity - in/hr. 

D = K/ d SM is the soil diffusivity - in/hr. , where 

d f 

d (SM )= moisture capacity i.e. the differential change in soil 

d | moisture with a change in capillary potential 



3 (SM) 

TT 



partial differential change in soil moisture 
with soil depth relative to the surface 
(when moisture supply is not a limiting 
condition) 



N.B. I 



saturated hydraulic conductivity K 



s at . 



AZS = AERATION ZONE STORAGE (UNSATURATED SOIL PROFILE) 



SM 


f 


y 

T 

p 


uns at . 


3 (SM) 


3 t 


3* 


§y 


SM 


K 
uns at. 



soil moisture - in. per unit depth 

capillary tension - in. 

gravity potential - in. 

total porosity of given soil profile 

unsaturated hydraulic conductivity. 
- Flow Equation - 

partial differential change in soil 
moisture with time, t 

partial differential change in soil 

moisture tension 

f (^) ; relates soil moisture as a function 

of the capillary tension ty and 

f (4»,SM) ; relates the unsaturated conductivity 

as a function of the capillary tension and soil 

moisture; these are secondary relationships 

required to be developed as aid in the numerical 

solution of the partial differential flow 

equation. 



IZS = INTERFLOW (INTERMEDIATE) ZONE STORAGE 



q 

Q c 
H 



V(KVH) 



szs 



sat. 
SATURATED ZONE STORAGE 



the outflow rate from the IZS - in/hr, , cfs 
relates to the volume of interflow in 
storage - cu.ft. 
total hydraulic potential. 

- Flow Equation - 

the Laplace's equation for unconfined steady- 
state flow in an isotropic porous medium 
saturated hydraulic conductivity. 



tf 



outflow rate from SZS - in/hr, , cfs 
volume in storage in base flow - in. 
total hydraulic potential (P/y + Z) , where 
P is the pressure potential, y the density 
of water and Z, gravity potential. 



SRO ■ SURFACE RUNOFF ROUTED INTO STREAMFLOW 



q - surface discharge {output) from the basin 
at time t, in cfs 



X = successive input volume up to time, t 
a 

h. ■ a time-area distribution parameter (distribution 

" a graph) for the successive input volumes - e, g, 

I.U.H. technique; t - summation over time t, 

I 
a = o 
for increments, a. 

INFL ■ INTERFLOW ROUTED OUT OF IZS INTO STREAMFLOW 

q^ ■ interflow rate - cfs at time, t 
t 

K = an interflow recession constant 

initial interflow rate (N. B, - routing 

relationship assumed to be similar to that 

of BFL). 



%% 



BFL m BASEFLOW ROUTED OUT OF SZS INTO STREAMFLOW 

q^ ■ base flow rate - cfs at time, t 
% 
X *■ a base flow recession constant - 

(N, B, - the routing function may require 

adjustments to account for variable 

component of ground -water recession) 

q = initial base flow rate, 
o 

TRO ■ TOTAL OUTFLOW AT BASIN OUTLET 

- Route via linear channel storage and flow translation 

I ss input 
Q - output 



K - storage coefficient 



dQ 

dt 



time rate of change in Q • 
ET ■ ACTUAL EVAPOTRANSPIRATION 



SSM - supply of moisture (water) available for 
evapotranspiration - amount available in 
IS, SS, IZS and AZS - in. 
R = radiation - langley. 
T = temperature - F, or C 
u = wind- m.p.h. 
H - relative humidity, H 
C-. = a crop (vegetation) -use factor. 
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Figure 4. HYDROLOGIC MODEL; A Schematization for an Integrated (Comprehensive) 

Hydrologic system Model 
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HYDRO LOGIC MODEL ; A Specific Model Approach 
A Simple Regression Model 

The following gives an outline for the definition of a 
parametric, linear regression model, according to the conceptual 
model scheme of Figure 5a, 

The basic relationship for a simple linear regression, 
rainfall-runoff model may be expressed as: 

R = a + bP • • • (1*3 

with b = Z (R. - B (P, " P)/ 2 (P. - P) 2 ... db) 

i-1 x x i-1 

and a - R - bP 

where R - runoff 

P = precipitation 

a = least-square regression constant (usually 

negative) 
b = least-square regression coefficient (usually 

< 1.0) 
1 « t, 2 r •..,<» 

Operation in the Whole System : One -Element Model, EQ - 

Figure 5a 

i) Relationship in terms of runoff, R: 
If P 1 C 



(2) 



o 
R - o 

If P > C 

I - % (P " C Q ) 

R = b P + a . - • ( 3 ) 

o o 

with 

a - -b C ■ constant* 



o 
'o 



where, C is a numerical parameter in EO, related to the surface 

o 

and subsurface storage or retention characteristic of the 
drainage basin, 

b is regression coefficient, an operator needed to convert P 
o 
(input to the system) to R (output from the system) and may be 
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related to the runoff characteristics of the basin; a is a regression 

constant, or a parameter related to the moisture deficit and 

retention of the soils (a minimum moisture requirement of the 

basin to be satisfied before the initiation of additional runoff 
by a given input P) . 

ii) Relationship in terms of basin loss, L: 



If P - C 



o 
L = P ... (4) 

If P > C 
o 

L , = P - R 

It m p - fa (P - C ) 

o o 

L = b' P + a* ... (5) 

o o 

with b* = (1 - b ) and 

o o 

a* - b C m constant, 
o o o 

b' is a regression coefficient or a parameter needed to 

convert portion of the P (retained or stored in the system) to 
loss L (output from the system) and may be related to the 
transpiration and evaporation characteristics of the vegetation 
and soil, respectively, in the basin, 

a' , is a regression constant, or a parameter related to 

a minimum moisture loss from the system when the moisture supply 
is not limiting. 

Operation in the Whole System: Three-Element Model, 
El, E2 and E3 - Figure 5a 

i) In Element El: 

If P - C ± 

R. = o ... (6a) 

L x = P ... (6b) 

If P > c 

1 

R- = P - C, . . . (7a) 
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where, C. is a numerical parameter in Element, El, related to 

the initial surface storage or retention characteristic of the 
system f 

L is the potential loss from Element, El, with a maximum 
capacity C. and 

R is the output from El which is a virtual input to Element, 
E2 

Element, El, may represent the surface storage characteristics of 
the basin (interception by vegetation is neglected) , 

ii) In Element E2: 

R = b-jR. • • • £i«J 

l 2 = (i - b x ) m x ■••« < 8b > 

substituting equation 7a into equation 8a and 8b: 

1 « fcL (P " C ± ) 

R m b P + & ± • •• ( 8c > 

with 

a 1 = -b C - constant, 

X 1 1 

L 2 = (1 - b x ) (P ~ C x ) 

L 2 m (1 " b l ) P + a 'l 

L„ - b\ P + a\ ••- (8d) 

with 

b< x = a - b x ) 

a ' = -b'_ CL = constant. 
1 11 

where, b, and b*. are factors, (regression coefficients) 
11 
needed for the conversion of input R, into runoff, R,or loss, 

L-, respectively. 

a 1 and a' are regression constants or parameters which may 

be related to the moisture deficit and moisture retention that 
has to be satisfied, before the initiation of runoff or loss 
from the system. 
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Element, E2 , may represent subsurface storage or retention 
characteristics of the basin. 

iii) In Element E3: 

L - L x + L 2 ... (9) 

where, L is the total potential loss from the drainage basin, 
i. e. that portion of the precipitation input that is not accounted 
for by the measured runoff. 

E3 is a 'dummy' element introduced to provide for the sum- 
mation of the potential losses from El and E2. 

Application of the Simple Regression Model 

Basic criteria ; 

a) - drainage basin; select a basin that exhibits a recurrent 

storage condition at the onset of each major seasonal 
change, e. g. onset of major spring runoff when the 
basin storage is at a maximum; or select a basin with 
runoff sequences that may be regarded as having 
negligible contribution from adjacent basins, 

b) - time interval; select a time interval for the analysis 

that indicates a negligible lag or autocorrelation 
effect in the respective time series of precipitation 
or runoff data, e.g. monthly, seasonal or annual series. 

Data Requirement 

- concurrent series of precipitation and runoff (monthly or annual 

series) , 
If monthly series are to be used, it is suggested that a 
normalization procedure be performed on the data prior to the 
regression analysis, e.g. use the logarithms of the original 
data assuming that the variates are log-normally distributed. 

Parameter Evaluation 

Gross estimates of C or C. may be assigned from knowledge 

of the major soil types in the basin and their area! distribution 
and moisture characteristics , the surficial geological charac- 
teristics, the recharge and discharge areas and the major aquifer 
storage characteristics. 

The constants and coefficients in the regression model 
equation (la) may be determined by operating with the concurrent 
series of P and R in a standard least-square regression analysis 
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subject to satisfying a threshold precipitation values 

R = a + bP, for P > P Q 

R = o, for P * P Q 

where, P is the threshold value {e.g. infiltration index) 
required to be satisfied before the initiation of runoff. 

Expanded (Multiple) Regression Model 

Consider the multiple regression: 

n 

y = I x k + a 

k = 1 K 

where, y is the dependent variable (runoff) 

X, are the independent variables (precipitation, evapotrans- 
piration, etc.); K = 1, 2.,.n independent variables and 
a is the regression constant. 

If the basin water balance or conceptual storage equation is used 
as a physical base , 



with 



where, R is the runoff and is a function of the precipitation, P, 
the basin loss, L, and the basin storage, S; in turn, L is a 
function of the evapo transpiration, E, and the loss to ground- 
water flow systems, G; and S is a function of the soil moisture 
storage, S , and the ground-water storage, S^. 

Operation in the Whole System: Expanded Three -Element 
(Multiple-Element) Model - Figure 5b 

i) Operation in Elements El, E2 and E3 

It should be noted that the operations in El, E2 and E3 
are as per previous development for the Three-Element Model, 
with the exception that the output R 2 from E2 represent a virtual 

surface runoff. Ihis is further augmented in E5 for base flow 
contribution to give output R from the system; 



R 


= 


f(P 


'• L* 


S) 


L 


= 


f l 


(E, 


G) and 


S 


= 


f 2 


(S a< 


V 



V - 6 

R 2 = b l R l *"° (10) 

where, the coefficient b, usually differs in value from the 

coefficient in equation (8a) of the Three-Element model. Also, 
the operations in E4a and E4b represent sub-operations extruded 
from El and E2 , 

ii) Operation in Element E4 

By definition, S = L, that is, input to storage is _ equivalent 
to the potential losses* accounting for evapo trans pi ration, change 
in storages and loss to regional flow system. 

Let C 7 be a numerical parameter in Element E4 related to the 
total water-holding capacity of the soil profile, i.e. total 
porosity, and assuming that the surface supply is an instantaneous 
input to title soil profile. 

If s <c 2 

S = S • « * (U a) 

S. = o • • • (lib) 

b 

If s >c 2 

S = S (capacity) ... (He) 

a a 

S, = S - S (capacity) 

where, S is the increase to available soil moisture in the soil 
a 

profile (unsaturated zone) and 

S, is the percolation (drainage) to the water table 

(saturated zone) . 

iii) Operation in Sub-Element, E4 a 

By definition S (capacity) = C.^, where, C 3 is a numerical 

parameter in E4a related to the moisture capacity of the 
unsaturated zone, after drainage by gravity (i.e. field capacity). 

a 3 

E ■ b«S . . « (12a) 

2 a 

AS = (1 - b~) S ... (12b) 

a d* a 

where, E is evapo trans pi ration (actual) 

feu is a parameter D r regression coefficient needed to 
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convert a portion of the available soil moisture to actual E, 

(l-b 2 ) is related to the soil moisture -re tent ion character- 
istics and 

AS is the incremental change in soil moisture storage. 

iv) Operation in Sub-Element, E4b 

Let C 4 be a numerical parameter in E4b related to a base 

storage condition of the local ground-water flow system, i.e. 
the aquifer storage condition necessary to sustain a minimum 
base flow (assuming that fixed proportions of the recorded 
precipitation are transformed simultaneously to base flow and 
changes in storage). 

If S K << C. 
D 4 

G = C. . . - (13a) 

4 

G = b 3 G -. • (13b) 

G = o (assumed) - • • (13c) 

AS, = (1 - b.) G ... (13d) 

D 3 

if s b * c 4 

G - S, + C. ... (14a) 

D 4 

G b = b 3 G 

G b = b 3 S b + a 3 * • • (14b) 

with a,, - b-C. - constant 

3 3 4 

G = o (assumed) 
and AS b = (1 - b 3 )G 

AS b = (1 - b 3 ) S b + a' 3 ... (14c) 

with a' ■ (1 - b 3 ) C 4 = constant. 

if s b »c 4 

G - S b + C 4 

G b = b 3 S b + a 3 
G r = b' 3 G 

G = b * 3 s b + a " 3 *•• (15a) 
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with a", = b'-C, - constant 

3 3 4 

and AS b = (1 - b 3 - b' 3 ) G 

AS b - (1 - b 3 - b» 3 ) S b + a'" 3 ... (15b) 

with a ,n 3 m {1 - b - b' 3 ) C 4 = constant, 

where, G is the total potential outflow from the local ground- 
water flow sy stern, 
G is the local ground -water contribution as base flow, 

G is the loss from the local ground-water flow system to the 
r regional flow system through deep percolation, 

b, is the regression coefficient needed to convert a portion of 

the ground -water storage to a base flow contribution and 
(1 - b 3 ) may be related to the aquifer storage characteristics. 

a., is a regression constant and may be related to minimum base flow 
contribution, 

b' is a regression coefficient needed to account for that portion 
of recharge lost, through deep percolation, to the regional 
flow system and 
(1 - b, - b' ) may also be related to the aquifer storage 

characteristics, 
a% may be related to a minimum deep percolation loss where 

recharge to local ground-water flow is not limiting, 
a' and a' " may also be related to the overall aquifer storage 

characteristics, and 
AS,, is the incremental change in ground-water storage. 

D 

v) Operation in Element, E5 

R = R 2 + G ... (16) 

where ,r is the total runoff from the system. E5 is a 'dummy' 
element introduced to provide for the summation of the conceptual 
surface runoff and base flow to produce the stream flow. 

Application of the Multiple Regression Model 

Consider the multiple regression model involving P, E and S 
as predictors of runoff, R: 

R - a + b,, P + b~E + b-S ... (16) 

o 1 2 3 

where, a , b 1 , P, b 2 , E, b 3 and S are as previously defined. 
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Definition of the Model Parameters 

Equation (16) can be written according to a regression model 

(Amorocho and Orlob, 1961) as equation (17) : 

m 

R = b, P. + b E + . I n a. P. + a^ ... (17) 

11 2 c 3 - I 3 1 o 

where, R is the observed runoff, 

P.. is the recorded precipitation, 

E is the computed estimate of evapotranspiration, 

Ea.P. represents the storage term, S, in which P. is 
3 J concurrent and antecedent precipitation 
values ; j = l,2,,,.,m period, 
a , a., b. , h 7 , are regression constants and coefficients. 

O lift. #■ 

The input parameters P and E are determined with a certain amount 
of error in their estimates viz,, 

a) Inflow - Index: 

P, = b'.P + e ... (18) 

lip 

where. P.. is an estimate of the areal average precipitation, 

P is the station (s) recorded amount, 

b' a station weight (coefficient) for areal distribution of P 

1 and 

e is the standard error of the average. 

P 

b) Evapotranspiration - Estimate: 

E = b' E ± e ... (19) 

c 2 c 

where, E is the areal average estimate of evapotranspiration, 

E c is the station (s) computed estimate, 

h % a station weight (coefficient) for areal distribution of 

E and e the standard error of estimate, 
c 

c) Storage - Index: 

If it is assumed that fixed proportions of the total recorded 
precipitation are transformed into ground-water outflow, sim- 
ultaneously producing changes in storage, then the storage index 
term of equation (17) may be interpreted as: 

AS + AS, + G = £*.?. + a ... (20a) 

a b r j j o 

where, AS , AS, are change in soil moisture and ground-water 
a storage and G is deep percolation loss. 

If G is assumed to be negligible, with an appropriate selection 
r of the time interval, the storage index, S, may be estimated 
with presumably small error, from: 
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S = I - G. ... (20b) 

where. 2L, is the net inflow at the soil surface and G. the 

base flow, 
which may be determined from total runoff hydrographs by hydrograph 

separation techniques, Incidentally, the net inflow may 

be examined on the basis of: 
m 

L T = E- a.P . + a + G. . . - (21) 

N j =1 j j o b 

Model Simplification 

The following statistical criteria are given for the sim- 
plification of the model relative to the selection of appropriate 
time periods for the analyses. At certain periods , some of the 
major parameters may tend to a constant due to recurring climatic 
conditions. For example, to examine the statistical degree of 
accuracy with which 1 can be estimated on the basis of a 
multiple regression model, 

n 
y - k | 1 x k + a . . . (22a) 

inhere, X, are random independent variables), 

usually the best choice of an end-point is when the variance of 

the observed values of a given variable is at a minimum,, or when 
its values are approximately equal to the average value for the 

period of observation. 

To examine this, consider the variance of the regressed 
estimate , y : 

Var (y) - E var (X. ) +2 I Co? (X.X. ) ... (22b) 

k K j,k 3 K 

where, X, and X, are specified independent variables. Minimization 

of the error in y may be examined by the following relationships: 

3 Var (y ) (23] 

I Var (X^) " ° "• U3) 

1 Var (y) = o ... (24) 

S Gov (XjX k ) 

that is, taking the partial derivative of the variance of y with 
respect to the variance of X, and the covariance of X, and X, , 

and equating to zero, respectively. 



(25) 
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As the sample estimate of Var (X k ) is given by: 

Var (X,) - 1 I (X 2 k - nX k 2 ) 
n - 1 k 

.*. IX 2 K = n X k 2 , for Var (X k ) - o ... (26) 

Hence X, -*• X, for minimum variance satisfying the partial 

derivative given by equation (2 3) . Similarly, the sample 
estimate of the covariance is given by: 

Gov (X.X.) = 1 1 (X.X, -nX.X.) ... (27) 



• • 



r^r k,j 

j"k 



ZX.X,. - n X.X k , for Gov OKA, J ■ o ... (28) 



Hence X.X, ■* X^ for minimum variance satisfying the partial 
derivative given by equation (24) . 

These criteria may be applied to model equation (16) . 

Consider a time interval such that the values of two of the 

independent variables approach their respective mean values, 
for each set of observations, 

Viz, E c - E c 

AS ■* AS 

or choose a period when the mean value of one variable becomes 
independent of the magnitude of the other variable, i.e. the 
conditional expectations of one variable* s values, given the 
other variable, tends to a constant: 

E (E | AS) ■* constant + If 

c c 

or 

E (AS|e ) + constant * AS 

A similar consideration may be given to other pairs of 
independent variables* 

By operating with the above criteria, the relationship 
(equation (16)) may be simplified: 

(a) for an observational period with a constant evapotrans- 
piration rate , 
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E +1 ^ C^ At; Var (E ) -* o ... (29) 

c e 1 c 

where, C_ is a constant evapotranspiration rate for At, the time 
interval (t, - t ) . 

1 Q 

(b) for an observational period with approximately equal changes 
in storage, - considering each storage component separately - 

i) for base flow: 

ASj+AS, * o; Var (AS ) ■* o 
and 6L-* G, = CLAt| Var (B^) ■* o ... (30) 

where, C, is a constant base flow rate for At, the time interval 

c% - t Q ) . 

ii) for a period with constant M (period when near maximum 

saturation or near complete depletion of moisture in the 
aeration zone can foe assumed) : 

AS + AS -* o ; Var (AS ) ■* o 
a a a 

Now considering together conditions (a) and (b) , the 
equation for surface runoff R may be expressed: 

R s = b l P l " (C E + C b } At "• (31) 

If E and G, cannot be estimated independently, then a 

general relationship between runoff and precipitation 
would result with regression coefficients that reflect 
the effect of evapotranspiration and base flow. 

In general, on a water-year basis with a carefully 
selected end-point, e. g, latter part of the drought 
season, when 

AS,-* AS, * o; Var (AS,)-* o 



and 



AS -*■ AS * o; Var (AS ) -*• o 
a a a 



G - o (assumed! 



the relationship may be expressed as: 

R - b, P- - E ... (32) 

11 c 
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If E is not estimated independently, the regression is: 
c 



R ■ b.P 1 + a x 



(33) 



in which case a regression would be the result with the effect 



of E reflected in a, and b, , 

C 1 A 



Data Requirement 

- historical record of the various variables. 

For an annual P-R model analysis (arranged on a yearly basis, 
beginning with each month) , a concurrent series is required! 



1 - P 

2 - R 

3 - E 



4 - 

5 - 



6 - P (snowpack storage) 



For monthly P-R model analysis, the following concurrent series 
is required: 



1 - P 

2 - R 

3 - E 

c 

4 - P (snow) 



5 - S. (at the beginning of the month) 

b 

6 - S (at the beginning of the month) 

a 

7 - S. (at the end of the month) 

b 

8 - S (at the end of the month) 

a 

9 - other meteorological parameters 

(temperature, humidity, radiation). 



Parameter Evaluation 



Estimates of G, , C« 



and C g may be assigned from knowledge of 



the major soil types in the basin and their areal distribution and 
hydraulic and hydrologic characteristics; estimates of C^ may be 

assigned from knowledge of the basin surficial geological char- 
acteristics , recharge and discharge areas and the major aquifer 
hydraulic and hydrologic characteristics, 

The constants and the coefficients of the regression models 
(e.g. equations 31, 32, 33) may be determined by least-square 
techniques through operations in the appropriate regression 
sub-routines. These regression sub-routines are the result of 
the appropriate application of the statistical criteria for the 
simplification of the multiple regression models equation 16 or 
17. 
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Whole System-, One-Element Model 



! 



EO 



Whole System: Three-Element Model 



r 



\ 



J El 

J 



h 



EO 



E2 H- 

/ I 



» 2 

E3 H 



\ 



J 



Symbols 

P - precipitation 
IL - net input 

R - runoff 

L, L , L 2 - loss 

EO,.,,E3 - element off 

the system. 



Figure 5a: HYDROLOGIC MODEL : A 'Definition Scheme fox a 

Parametric, Linear Regression Model 
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Whole System: Expanded Three-Element 

( Multiple Element) Model 



E 



r 



LJ 



El 



\ / 

\ s 



E3 L 



T 

S = 



• N 

/ \ 

I E2 . h 



* E4a v »« a / E4 1 

v J I I 






EO 



/ x 

' E5 L 

I / 



« s b f VAU \ % 

__*J E4b j 



"1 



R 






Symbols 

P - precipitation 
IL - net input 

R,~,R- runoff 
2, 



1 

S 

S 
a 

s b 
% 

G 



L , L - loss term 

- storage 

- soil moisture storage increment 

- ground-water storage increment 

- base flow contribution 

- loss to regional flow 



E - evapotranspiration loss 
E1,...,E5 - element of the system 



Figure 5b: HYDROLOGIC MODEL: 



A Definition Scheme for an Expanded, 
Parametric, Linear Regression Model 
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Figure 5c : HYDROLOGIC MODEL: Flow Chart for Parametric* 

Linear Regression Models 



Whole System- One-Element Model, EO 
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' 
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INPUT DATA: 
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PARAM,. , EQN, - 
Coef. fi Const. 

; r a ,a* f b ,b" 

o o o o o 


















1 


} 






- 


It 

R 

I 

i = 1, N 
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no y^ 


< C 

- O <s 


s^ ye- 


I 


R. = o 






\ 1 - 




L i = p i 


















— ►> 


R. = b (P. - C ) 
i o i o 

L. = P. - R. 












D. = R. - R. 

ISE= tt(R.- Rjj 2 } 
{ 1 j 








OT 












R. ■ b P. + a 
1 ox o 

L. = b' P. + a" 
x o i o 




1- - w 










1 


' 












OUTPUT 

PRINT: P ± , R ±f 
R , T. . _~n ... TSE 










STOP 
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Whole System: 'Three-Element Model, E1,E2 & 13 



INPUT DATA: 
P. R. 



START 



Si 



PARAM, , EQN, - 
Coef. & Const, 

C l' a l' a, l'' b l' b 'l 



El 



E2 







no 


^ 


■^P. < c, ^x yes 




! 


r 






• 


' 






R li " V C l 




R 1 . -o 

L li ■ P i 
Lli - 










__ . _. _ - _ 




— • 




— — 




. 





— 


— *■ 


R. = b 1 H li 








L i =L li +L 2i 




E3 












1 










or 








\ 


? 








— * 


R. =b lF . ♦ a x 
L 2i = b 'l P i + ' 


''l 






D. = R. - R. 




OUTPUT 
PRINT: P, ,R, ,R. 

L li' L 2i' L i 

D. , ISE 










ISE 


■ ( MR i- B i ) } 










1 vfc 2 1 












11 "x 
















! 


' 






STOP 
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Whole System 



Expanded Three-Element (Multiple-Element) 
Model, El, E2, E3 , E4 , (E4a, E4b) , E5 
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INPUT 
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S, - L, 
l l 



no 




bi l 3 



S , 9 S, 
ax x 

S, . = o 
bi 



E. = b„S . 
1 2 ai 



AS .«(l-b_)S 
ai 2 ai 




G. = S. . + C. 
l bi 4 






G, , = b_G, 
bi 3 x 






G . - o 
rx 






As^.^X-bjG. 
bx 3 x 






or 






G. . = b S + a 
Tsi 3 bi 3 






AS. .= (l-b.)S, . + a* 
bx 3 bx 


3 





G, = S. . + C. 

x bx 4 

G, .« b,G. 

bx 3 x 

G .= b'G. 
rx 3 x 

As. . = (l-b.-b')G. 
bx 3 3 x 



G, . = b_S. . + a, 
bx 3 bx 3 

G . - b'S, . + * rt , 
rx 3 bx 3 

As,.= (l-b.-b' )5, .+ 
bx 3 3 bx 



_£££_ 



G. = C 
x 4 

G . = o 
rx 

GL . = b,G. 

bx 3 x 

As. .*(l-b_)G. 

bx 3 x 
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ise= CEC* A -ft^ } h 

1 Ei.2 1 } 



OUTPUT 

PRINT: P. f JL ,R. ,L. 

l i ll 

D. , X5E 

S i' S ai' S bi 

V AS ai' AS bi 

G i ,G bi' G ri 



STOP 
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LIST OF SYMBOLS IN MODEL (FLOW CHART ) 

ONE- ELEMENT MODEL BO 

P, ■ recorded precipitation - in, 

R, = recorded runoff - in, 

C = parameter for basin retention - in. 
o 

a . a' = regression constants 
o' o y 

b , b' = regression coefficients 
o* o " 

i = 1, 2, 3,....N periods of observation 
It. - computed runoff - in. 
L. = computed loss - in. 

D. - regression residual - in. 
ISE = integral square error 

THREE -ELEMENT MODEL El, E2, & E3_ 

P, = recorded precipitation - in, 

R. - recorded runoff - in. 

C, «■ parameter for basin retention - in. 

a, , a' 1 - regression constants 

b 1 r b 1 , - regression coefficients 

L,. « a loss term, from El - in. 

R 1 , ■ computed net input to E2 - in, 

L-. = a loss term, from E2 - in, 

L. ■ sum of the losses - in. 

R. = computed runoff -in, 

i = 1, 2, 3....N periods of observation 

MULTIPLE-ELEMENT (Expanded 3-Element) MODEL 

Element E4 , E4a, E4b & E5 (see definitions above plus) 

S. = storage input (= potential total losses) - in. 

C- = a parameter for basin retention - in. 

CL * a parameter for basin retention - in, 

S . - soil moisture storage increment - in, 

ai 
S . = ground-water storage increment - in, 

E. ■ evapotranspiration - in. 
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AS . = change in soil moisture storage - in. 

C- - a parameter for basin storage - in, 

6. = a potential outflow increment - in, 

G. . = a baseflow contribution - in, 

G . « loss to regional flow system - in. 

AS. j = change in baseflow storage - in, 

R„ . - a surface runoff component - in. 

R. = computed runoff - in. 

a.,, a 'v a "v a " l = re 9 ress i° n constants 

b~# b_, b' = regression coefficients 
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HYDRO LOGIC MODEL : 

A Brief on Stochastic Models 

The behaviour of a hydrologic system and/or process varies 
with a sequential time function. This variable process has a 
degree of uncertainty which can be described by probability laws 
and, as such, the system or process may be regarded as stochastic 
(an index family of random variables) . That is, there is an element 
of chance in the values which the variables may take in the 
interrelated hydrological processes which define the average 
behaviour of the system. 

Parametric hydrology generally deals with discrete hydrologic 
events, such as storms or floods. A parametric approach is 
deterministic in the sense that derived functional relationships 
between the input-output processes attempt to define the behaviour 
of the system for all times; that is, a given set of causes always 
produces the same set of effects. Stochastic hydrology, on the 
other hand, concerns itself with the time -sequential properties of 
the discrete events. It attempts to formulate a mathematical 
model for an observable phenomena which changes in time, in a way 
that is not completely predictable, but is expressible in terms 
of probabilities. 

A conventional approach to stochastic hydrologic system 
modelling usually takes the form of a hybrid model (deterministic 
and stochastic). To achieve this hybrid model, a firm foundation 
based on the deterministic model concepts is required. The model 
may be used to transform the recorded series input of a variable 
with known or assumed probability distribution to values of an 
output with unknown probability distribution. Time-series analysis 
may then be applied to examine the generated output for its 
stochastic properties. Alternatively, the stochastic properties 
of the input variable may first be determined, then a deterministic 
relationship is used to transform the input to the required 
output. The generated output is presumed to have stochastic 
properties similar to that of the input variable. 

Another alternative method, such as in the case where a 
deterministic relationship is known to exist between two time 
series, but cannot be specified explicitly, is to establish 
through regression or correlation, empirical relationships between 
parameters derived from statistical analyses of the individual 
time series. These empirical relationships are then used as a 
deterministic base, for the transformation procedures. 

A more aggressive approach (Chow & Kareliotis) to stochastic 
hydrologic system modelling, Figure 6, is to formulate a stochastic 
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model for the hydrologic system in which runoff is considered as 
the integral of three component stochastic processes, viz., 
change in conceptual watershed storage, the total rainfall input 
and the total watershed losses, mainly evapotranspi ration. In 
the analysis, the integrated stochastic processes are treated 
as a three-dimensional vector or multiple-time series with each 
time series consisting of a deterministic and a random component, 

The mathematical notation given for each stochastic process 
in Figure 6 is called a time series. For example the set of 
observations (P (t) | tET) is called a time series for precip- 
itation; i.e. for each t in T the observation P (t) is an 
observation of the random variable. Prior to analyzing the time 
series, it is necessary to .assume a model for the stochastic 
p ro ce s s , e , g . mo v ing ave ra ge , MA , ha rmo n ic , H » a uto - re gres s io n , 
AR, or Markov Chain, MC. Standard statistical techniques in 
time-series analyses (correlogram and spectral density analysis) 
add in the identification of the model inherent in the data 
sequence,, through goodness of fit tests on the assumed 'model by 
comparing generated output with historical data sequences. 



* 






Precipitation - Simulator 

- available data - historical, discrete time 
series : 

* {p (t)i t et! 

- assume that fe (t) j t £T> is a stochastic 
process 

- assume a type of model (moving average, 
harmonic, or Markov Chain, etc.) 

- analyse time series (correlogram or spectral 
density analysis) 

- identify and choose appropriate simulator 
(Model). 
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Evapotransplration - Simulator 

- available data - computed estimates, 
discrete time series: 

f E (t) ; t £T} 

- assume that' {e (t) j t €T> is a 
stochastic process 

- assume type of model 

- analyse time series 

- identify and choose appropriate 
simulator. 



Notations : 

P (t> - Precipitation observation 

E (t) - Evapotranspiration observation 

S (t) - Conceptual watershed storage 

R (t) - Total runoff observation 

t ET - Denotes the discrete time t, 

belongs to (sub-set) of an 

interval of real time T? 

T - o < t £ I, 

where L is the length of the 

time interval. 



Concept ual Watershed Storage - Simulator 



- available data - computed estimates, 
discrete time series: 

fe (t) $ t ET} 

- assume that fe {t)j t ET} is a 
stochastic process 

- assume type of model 

- analyse time series 

- identify and choose appropriate 
simulator. 



Transformation System 

- assume multiple time series: 
(P(t), E(t), S(t), R(t) ; t ET} 

- assume that the multiple time series 
is a stochastic process 

- assume type of model (mass balance) : 
S(t + 1)- s(t)= P(t) - R(t) - Eft) or 
E(t)= P(t) - R(t)-(S(t + 1) - S(t>) or 
R(t)« P(t) - E(t)-(S(t + 1) - S(t)) 

- optimize model parameters. 



Runoff - Simulator 

- available data - historical, 
discrete time series: 

(R 111 ; t ET } 

- assume that (r (t)i t ET> is a 
stochastic process 

- assume type of model 

- analyse time series 

- identify and choose appropriate 
simulator, 



Figure 6. HYDBDLOGIC MODEL t 



A Schematization for a Type of Stochastic 
Hydro logic System Model. 
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